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Abstract 

We discuss the recently established equivalence between the Laplacian growth 
in the limit of zero surface tension and the universal Whitham hierarchy known 
in soliton theory. This equivalence allows one to distinguish a class of exact so- 
lutions to the Laplacian growth problem in the multiply-connected case. These 
solutions corerespond to finite-dimensional reductions of the Whitham hierarchy 
representable as equations of hydrodynamic type which are solvable by means of 
the generalized hodograph method. 
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1 Introduction 



Dynamics of a moving front (an interface) between two distinct phases driven by a har- 
monic scalar field often appears in different physical and mathematical contexts and has 
a number of important practical applications. Processes of this type are unified by the 
name of Laplacian growth. Their key common feature is a harmonic field which serves 
as a potential for the growth velocity field. 

The most extensively studied is the case of two-dimensional spatial geometry. To be 
definite, we shall speak about an interface between two immiscable and incompressible 
fluids with very different viscosities on the plane. In practice the 2D geometry is realized 
in a narrow gap between two parallel glass plates (the Hele-Shaw cell). In this version, 
the problem is also known as the Saffman- Taylor problem. For a review, see [1]. 

To be more precise, consider a compact plane domain occupied by a fluid with a 
negligible viscosity (water). We call it water droplet. Let its exterior be occupied by a 
viscous fluid (oil). The oil/water interface is assumed to be a smooth closed contour. 
Oil is sucked out with a constant rate through a sink (a pump) placed at some fixed 
point (which may be at infinity) while water is injected into the water droplet. More 
generally, one may consider several oil pumps and several disconnected water droplets. 
In this general setting the problem was discussed in [2] -[6]. 

In the oil domain, the local velocity V = (V x , V y ) of the fluid is proportional to the 
gradient of pressure p = p(x, y) (Darcy's law): 

V = -nVp 

where k is called the filtration coefficient. It is inversely proportional to viscosity of the 
fluid. For future convenience, we choose units in such a way that k for oil be equal to 
j. In particular, the Darcy law holds on the outer side of the interface thus governing 
its dynamics: V n = — \d n p. Here V n is the growth velocity, which is by definition normal 
to the interface, and d n is the normal derivative. Since the fluids are incompressible 
(Vy = 0) the Darcy law implies that pressure p is a harmonic function in the exterior 
(oil) domain except at the point where the oil pump is placed. There is a logarithmic 
singularity at this point. 

So, the pressure field is a solution of the time-dependent boundary value problem for 
the Laplace operator A = d 2 x + dy with a prescribed logarithmic singularity and with 
certain condition on the boundary of the growing domain. This condition is fixed by 
the following reasoning. If viscosity of water is small enough, pressure is uniform inside 
water droplets (p = const). Hereafter, we assume that pressure does not jump across 
the boundary, so the function p(x,y) is constant along it. This means that the surface 
tension is set to be zero. To summarize, the mathematical formulation of the Laplacian 
growth, in the limit of zero surface tension, is: 

V n = ~\d n p, 

Ap = in oil with dp = along the boundary, (1) 
p(x,y) = logp 2 + ... as p — > 

Here p is the distance between the sink and the point x, y. If the sink is at infinity, 
what is often assumed in the outer growth problem, the last condition is replaced by the 
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asymptotics p = — log(x 2 + y 2 ) + . . . very far away from the droplet. If the oil domain is 
multiply-connected (i.e., there are more than one water droplet), this formulation should 
be supplemented by some additional physical conditions for pressure in the water droplets, 
which we do not discuss here. The problem consists in finding the time evolution of the 
interface between the fluids. 

The simply-connected case (one water droplet) allows for an effective application of 
the conformal mapping technique (see, e.g., [1, 5]). Passing to the complex coordinates 
z = x + iy, z = x — iy, one may describe the growth process in terms of a time dependent 
conformal map z = f(w,T) from a fixed domain of a simple form, say the exterior of 
the unit disk, onto the evolving oil domain. The interface itself is the image of the 
unit circle w = e 1 ^, < <fi < 2tt. The dynamics (1) is then translated to a nonlinear 
partial differential equation for the function f(w,T), referred to as the Laplacian growth 
equation. If the oil sink is at infinity, it reads 



It first appeared in 1945 in the works [7] on the mathematical theory of oil production. 

Later, some exact finite dimensional reductions of the Laplacian growth equation 
were found. Asuming that the derivative of the conformal map / is a rational function, 
it reduces, upon a direct substitution, to a finite dynamical system for parameters of this 
function (say, for its poles, zeros or residues). See [8], [9] for different versions of such 
dynamical systems. Moreover, it turnes out that these systems can be actually integrated 
in the sense that all time derivatives are removable. As a result, the poles are represented 
as implicit functions of time. 

Existence of finite-dimensional reductions was a serious motivation for searching an 
underlying integrable structure of the model. In case of one water droplet, such a struc- 
ture was identified in [10] to be the 2D Toda hierarchy in the limit of zero dispersion. The 
Laplacian growth equation is the string equation of the Toda hierarchy. The conformal 
map plays the role of the Lax function. 

The vector field d/dT in the space of domains defines a particular flow of the Toda 
hierarchy. Other flows are "frozen" until oil is sucked by a single pump only. To "unfroze" 
them, one should allow oil pumps to work at any point in the plane. In general there are 
as many flows as types and locations of oil pumps. It is convenient to associate to each 
pump its own time variable which is the total amount of oil sucked by the pump. It is 
known [6] that the result of the Laplacian growth evolution, in the case of several pumps, 
depends only on the total amount of oil sucked by each of them and does not depend on 
the order in which they operate. In the context of integrable systems, this fact reflects 
commutativity of flows of the Toda hierarchy. 

Recently, this approach was extended [11] to the Laplacian growth of several water 
droplets. In this case the dynamics turns out to be equivalent to the universal Whitham 
hierarchy on the moduli space of genus g Riemann surfaces, where g is the number of 
water droplets minus 1. In full generality, this hierarchy was introduced by Krichever 
[12, 13] in the context of slow modulations of exact periodic solutions to soliton equations. 
The dispersionless Toda hierarchy is its simplest particular case at g — 0. The universal 
Whitham hierarchy is known to have a lot of meaningful finite-dimensional reductions 
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Figure 1: Water droplets D a (g = 3). 

for any g. They are systems of differential equations of hydrodynamic type which can 
be implicitly solved by the generalized hodograph method [14]. Applying these results 
to the Laplacian growth problem, we obtain a large family of solutions in the multiply- 
connected case. The infinite dimensional contour dynamics for them is reduced to finite 
systems of differential equations which can be integrated. 

The aim of the present paper is to study this class of finite-dimensional reductions 
and to clarify their algebro-geometric meaning. 

2 Laplacian flows and the universal Whitham hier- 
archy 

In this section we explain the precise meaning of the statement that the Laplacian growth 
is described by the Whitham equations. 

We use the following notation (Fig. 1). Let D oil be the region of the plane occupied by 
oil (a non-compact domain containing infinity), and D be the region occupied by water. 
We assume that there are g + 1 water droplets, which are compact domains bounded by 
smooth non-intersecting closed contours 7 a , a = 0, 1, ... ,g. Let D a be the a-th water 
droplet, so that D is their union. By 7 we denote the union of the 7 Q 's: 7 = U^, =0 7 a . 
Our assumptions imply that no contour 7 Q lies inside another. 

2.1 Commutativity of Laplacian flows 

Any growth process defines a flow in the space of contours 1 , or, in the infinitesimal 
version, a vector field. Integral curves of this vector field are just one-parametric families 
of contours 7(T) obtained from any initial contour 7(0) by the evolution in the time T. 

Flows corresponding to the Laplacian growth (Laplacian flows) are determined by 
location of oil pumps. Namely, given an oil sink at a point a e D n, the solution of the 
boundary value problem 

Ap(z) =4tt5 (2) ( z — a) , dp = on 7 (3) 

1 Depending on a particular problem, the flow can be only defined on a proper subvariety or a sub- 
domain in the space of contours, but this is not important for us here. 
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in D oil (here 5^ 2 \z — a) is the two-dimensional delta- function supported at the point a), 
with some additional conditions on the boundaries of D Q 's in case of more than one water 
droplet (see [11] for details), uniquely defines an infinitesimal deformation of 7 by setting 
the vector of normal velocity to be — \d n p. Such a deformation is a tangent vector to the 
space of contours. For brevity, in this paper we do not discuss other types of Laplacian 
flows considered in [11]. 

Given a Laplacian flow, a natural parameter along the trajectories in the space of 
contours is the amount of oil, T, sucked by the pump which defines the flow. Since the 
both fluids are assumed to be incompressible, T is simply the change of the total area of 
the water droplets. (If the pumping rate is constant, T can be also identified with time.) 
Let dr be the corresponding vector field. It is defined in an invariant way, i.e., it does 
not depend on the choice of local coordinates in the space of contours. 

As is common, we identify vector fields with derivatives along them. The derivative 
of any functional F on the space of contours along the vector field (the Lie derivative) is 
defined as follows: 

a r F = . im F ^ r + £ »-™ 

If there are several sinks (at the points a, b, etc), we write d T ( a ), <9 T o), etc for the 
corresponding vector fields. 

The notation is justified by the following theorem [6]. 

Theorem 1. The vector fields d T ( a ), d T (b) commute for all a,b e D oil . 

It is an infinitesimal formulation of the fact that the result of the Laplacian growth 
evolution (i.e., the shape of D oil ) depends only on the amounts of oil T^ a \ sucked 
from the points a, b and does not depend on other details of the process (assuming that no 
topological change or singularity occurs). In its turn, this fact follows from the reduction 
of the Laplacian growth to the inverse potential problem in two dimensions and from 
local uniqueness of solutions to the latter. For more details see [6]. 

This result suggests to treat d T ( a ) as a partial derivative in a space with coordinates 
T( a \ To construct a simple example of such a space, fix N points aj and consider the 
variety of contours which can be obtained from some initial configuration as a result of all 
Laplacian growth processes with oil pumps at the points cij. The resulting shape of the 
droplets is uniquely determined (if no singularity occurs) by amounts of oil sucked out 
of each point a,j, i.e., by the quantities This configuration space is N- dimensional 

and are local coordinates in it. 

2.2 The Schottky double 

Our next task is to carry the Laplacian flows over to compact Riemann surfaces with a 
real structure. To do that we need the Schottky double construction [15] which we recall 
here. 

With the domain D n (as well as with any plane domain with a boundary) one canon- 
ically associates its Schottky double 

E = D oi i U 7 U D oil 
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where D* is another copy of D oil (its "opposite side" ) attached to it along the boundary 
7. In what follows, we shall call D oi i and D* n the front side and the back side of the 
double, respectively Formally speaking, a point P G E is a pair P = (z,o~), where 
z G D oi i and o = ± indicates front or back side. If z G 7, the pairs (z, ±) represent the 
same point. There is a natural projection E — > D oil which sends a point P — (z,a) & E 
to the point 2 G D oil . In the sequel, we sometimes write simply z for (z, +) and z* for 

(*,")■ 

By adding a point at infinity for each copy of D n, E becomes a compact surface 
naturally endowed with a complex structure. The complex structure on the front side 
is the same as in D n while the back side has the opposite complex structure. In other 
words, the holomorphic coordinate is z on D oi i and z on D* n . Whence E is a compact 
Riemann surface with the antiholomorphic involution P 1— > P* which sends (z, a) to 
(z, —a). The two sides of the double are in a sense mirror images of each other, so we 
shall call points connected by the involution mirror-symmetric points. Boundaries of the 
droplets are fixed points of the involution. A schematic view of the Schottky double see 
in Fig. 2 below. 

If the contours 7 a are smooth, the double is, by construction, a smooth Riemann 
surface of genus g. Indeed, each point of E, including points of the boundary of D n, has 
a disk-like neighborhood. 

Note that the set of fixed points of the involution consists of exactly g + 1 closed 
contours, which is the maximal possible number for surfaces of genus g. Riemann surfaces 
with antiholomorphic involution having this property are called real. 

In what follows it is implied that a canonical basis of a and b-cycles on the double 
is fixed. This can be done in many ways. For example, one may take a a -cycles to be 
boundaries of the droplets D a for a = l,...,g and take b a -cycles to be some paths 
connecting the boundary of the 0-th droplet with the a-th one on the front side and 
going back through the opposite side of the double. Each choice of the basis of cycles 
corresponds to certain physical conditions in the water droplets which are discussed in 
detail in [11]. For our purpose here it is enough to imply that one or another basis is 
fixed. 

Any pair of meromorphic functions </?, (p (which are functions on D n and on the 
complex congugate domain respectively) such that <p(z) = <p(z) on the boundary defines 
a meromorphic function on the double. Namely, this function takes the value <p(z) at 
the point z on the front side and the value <p(z) at the mirror-symmetric point z* on the 
back side. Similarly, a meromorphic differential on the double is represented by a pair 
of meromorphic differentials <p(z)dz and (p{z)dz such that ip(z)dz = <p(z)dz along the 
boundary contours. 

If a meromorphic differential du(z) on D n is purely imaginary along the boundary, 
the analytic continuation to the opposite side of the double is especially simple. Indeed, 
in accord with the Schwarz reflection principle, its meromorphic extension to the back 
side is — du)(z), so for each pole of such a globally defined differential there is a mirror 
pole on the opposite side. In this case we shall say that dw is a symmetric differential on 
E. 
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2.3 Abelian differentials on the double 

To each pump we assign an Abelian differential on E. On the front side it is the differen- 
tial d z pdz, where p is the pressure field generated by the pump in the viscous fluid. Since 
dp = d z pdz + d z pdz = on the boundary, d z pdz is purely imaginary along the boundary 
So, its analytic continuation to the back side of the double is —d z pdz. Conversely, given 
a symmetric Abelian differential on E, we place pumps at its poles on the front side. 
Strengths of the pumps are determined by the residues. 

Note also that any symmetric Abelian differential dQ on E defines a tangent vector 
to the space of contours via setting the normal velocity at any point of the contour to be 

V n = - l -d n Re fdn 

If poles of dfl are fixed in the z-plane, this defines a vector field corresponding to a 
Laplacian flow. 

For example, if oil is sucked (with a unit rate) from a point a G D oi i, then 

p(z) = 2G(a,z) 

where G is the Green function of the Dirichlet problem in D ii. The Green function 
is defined as a solution of the Laplace equation in D oil with the logarithmic singularity 
of the form G(a,z) — > log\z — a\ as z — > a and such that dG = on the boundaries. 
Therefore, the differential dW ( - a \z) = d z pdz has a simple pole at z = a with the residue 
f and no other singularities in D Q n. The analytic continuation of this differential to the 
back side is — dW^ a \z) which gives a simple pole at the mirror point a* with the residue 
— I. The differential constructed is a symmetric Abelian differential of the third kind (a 
dipole differential): 

!2d z G(a, z)dz on the front side 
(4) 
—2d z G(a, z)dz on the back side 

Remark. For g > 1 these conditions do not yet define the differential uniquely. To 
fix it, a normalization is needed. We normalize the differential dW^ by the condition 
that it has zero a-periods. Depending on the choice of cycles, this means certain physical 
conditions: for instance, fixed injecting rates or fixed pressures in the water droplets. See 
[If] for details. 

2.4 The Whitham hierarchy 

The Whitham hierarchy structure reveals itself in evolution of the Abelian differentials 
under Laplacian flows E — > E{T). 

The holomorphic coordinate z on the front side of E is independent of time by con- 
struction. It defines a connection in the space of real Riemann surfaces of genus g. 
Correspondingly, the Abelian differentials become time dependent: dW = dW(z,T). At 
any fixed z, dW(z) is a functional on the space of contours, so one can take derivatives 
along Laplacian vector fields. 
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In [11], it was shown that the Laplacian growth of D oil , with zero surface tension, 
is equivalent to the system of Whitham equations for symmetric Abelian differentials of 
the third kind dW^ on the Schottky double of D oil . The main result is the following 
theorem [11]: 

Theorem 2. For any pair of pumps (at the points a,b G D oily ) 7 the following relations 
hold: 

d T(a) dW {b) = d T(b) dW {a) (5) 

This statement follows from (4) and the Hadamard variational formula [16] for the Green 
function G(z,a). The Hadamard formula implies that d T ( a )G(b,c) is symmetric under 
permutations of the points a, b, c. 

The theorem allows one to identify Laplacian flows with Whitham flows on the (ex- 
tended) moduli space of real Riemann surfaces. Equations (5) constitute the universal 
Whitham hierarchy [13]. In the soliton theory, they are obtained by averaging solutions 
of soliton equations over fast oscillations [17]. 

Since the Laplacian vector fields commute, one may look for a general solution to 
equations (5) in the form 

dW {a) = -d T( a)dtt (6) 
where dQ is called the generating differential. 

Remark. In general, singularities of the generating differential may be of a very com- 
plicated nature. Equations (6) mean, however, that all of them do not move under the 
Laplacian flows. In other words, they are integrals of motion. The only parameter that 
changes in the course of the Laplacian growth with a sink at the point a is the residue of 
dQ at the pole at a. This means, in particular, that the Laplacian flows with point-like 
pumps can only add simple poles to singularities of dQ. 

In the next section we consider an important special case when all singularities of the 
generating differential are poles, i.e., dQ is meromorphic on E. 

3 Algebraic domains 

Families of special solutions to the Whitham equations with a meromorphic generating 
differential were constructed in [12]. In [11], it has been shown that they describe Lapla- 
cian growth of algebraic domains. It is clear that an algebraic domain remains algebraic 
under Laplacian flows with point-like pumps. This fact allows one to reduce the infinite 
dimensional configuration space of contours to a finite dimensional one. 

3.1 The Schwarz function 

For a more explicit characterization of the class of algebraic domains that we are going 
to deal with, consider the analytic continuation of the function z defined in D oil to the 
back side of the Shottky double. 
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Q 



E 




Figure 2: Poles of the functions Q and E on the Schottky double (g = I). The upper 
and lower halves of the torus give a schematic view of the front and back sides of the 
double, respectively. 



Definition. The domain D oil is said to be algebraic if the function z is extendable to a 
meromorphic function on the Schottky double o/D u. 

To clarify the meaning of this definiton, we recall the notion of the Schwarz function [18]. 
Given a bounded plane domain, the Schwarz function is the analytic continuation of the 
function z away from the boundary. Note that for such an analytic continuation to exist, 
the boundary needs to be analytic, not only smooth. By definition, the Schwarz function, 
S(z), is a function analytic in some neighborhood of the boundary contour (contours) 
such that 

S(z) — z on the contour 
If it exists, the Schwarz function must obey the "unitarity condition" 



S(S(z)) = z 



which means that the function S(z) is inverse to S(z). In general, the analytic continu- 
ations from different components of the boundary of a multiply-connected domain may 
lead to different functions. Algebraic domains correspond to a rather special situation 
when they lead to a common function which is well-defined everywhere in D n. 

If the Schwarz function exists, the analytic continuation of the function z to the 
back side is just S(z). Therefore, a plane domain is algebraic if and only if the Schwarz 
function is extendable to a meromorphic function in it. 

Let P G E be a point on the Schottky double of an algebraic domain. Consider the 
function Q(P) on E defined as follows: 



Q(P) 



S(z) if P is on the front side 
z if P is on the back side 



Here z = z(P) is the ^-coordinate of the point P = (z, a) on E. By definition of 
the Schwarz function, Q(P) extends continuously across the contour 7 and so it is a 
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well-defined meromorphic function on E. Along with the Q(P), we introduce a "mirror- 
symmetric" function 



E{P) = 



if P is on the front side 



S(z) if P is on the back side 

which is meromorphic on E as well. It provides the above mentioned meromorphic 
extension of the function z to the whole double. Note that the functions Q, E are 
connected by the relation 

S(E(P)) = Q{P) (7) 

and transform one into another (with complex conjugation) under the antiholomorphic 
involution: 

Q{P*) = E{P) (8) 

Suppose that the Schwarz function in D oi i has poles at some points pj G D oi i with 
multiplicities rij. Hereafter, the points pj are assumed to be in general position. Set 
N = J2j n j- Then the number of droplets can not exceed N + 1, or, what is the same, 
genus of the Schottky double is less or equal to iV: 



9<N = J2 



n 



This rough estimate 2 follows from the fact that any nontrivial function on a Riemann 
surface of genus g has at least g + 1 poles counted with multiplicities (if everything is in 
general position). Applying this statement to the function Q, we get the estimate. 

To summarize, on the Schottky double of an algebraic domain there exist two mero- 
morphic functions connected by the involution as in (8). The notation Q, E is chosen to 
emphasize the similarity with the approach developed in [13]. 

Remark. The class of algebraic domains introduced here and in [6] coincides with 
the class of "quadrature domains" from [19]. The class of domains corresponding to 
"algebraic orbits" of the Whitham equations [13] is broader. It includes domains for 
which the generating differential is meromorphic on E with "jumps" across certain cycles. 



3.2 A complex curve associated with an algebraic domain 

For algebraic domains, the Schwarz function is an algebraic function, i.e., it satisfies 
a polynomial equation R(S(z), z) = 0. The easiest way to see this is to consider the 
analytical continuations of the functions z, S(z) to the whole double, i.e., the functions 
E,Q. 

If S(z) has N poles in D oi i (counting multiplicities), then the function Q has N + 1 
poles on E (N poles on the front side and a simple pole at infinity on the back side). 
Similarly, the function E also has N + 1 poles. Therefore, we have two meromorphic 

2 This estimate is rather rough indeed. It can be improved by taking into account the reality condi- 
tions. For example, in the case N = 1 our inequality gives g < 1 while a more detailed analysis (see, 
e.g., Theorem 7 and Corollary 7.1 from [19]) shows that g = 0. 
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functions with N + 1 poles on a smooth Riemann surface. A simple theorem of algebraic 
geometry states that these functions are connected by a polynomial equation of degree 
N + 1 in each variable: 

R(Q,E)= ]T A nm Q n E m = (9) 

n,m=0 

Symmetry under the involution implies that the matrix of coefficients is Hermitian: 
A nm = A mn . Restricting the relation between Q and E to the front side of the dou- 
ble, we get the polynomial equation 

R(S(z),z) = 

The unitarity condition for the Schwarz function is equivalent to the hermiticity of coef- 
ficients of this polynomial. 

The equation R(z, z) — defines a complex algebraic curve r which may serve as 
another realization of the Riemann surface E. The boundary contour is recovered as the 
set of real solutions x,y to the equation R(x — iy, x + iy) = 0. It is the real section of 
the curve r. 

Remark. We see that the boundary of any algebraic domain is always described by a 
polynomial equation. The converse is not true. An arbitrary polynomial equation with 
hermitian coefficients in general does not define an algebraic domain. 

Despite of the fact that E is a smooth surface, the curve r is not smooth for N > 0. 
It contains singular points (double points in case of general position). This is seen from 
comparing genus of the double with the upper bound of genus of the curve r. For a 
generic smooth curve of the form (9) genus equals N 2 , so the fact that the genus does not 
actually exceed N just gives evidence for the presence of singular points on the curve r. 

How do the singular points emerge? Consider the map P i— > (Q(P), E(P)) from E 
to r which takes a point P to the pair of complex numbers (Q(P), E(P)) connected by 
the polynomial equation. If two different points of E are mapped to one point of r (i.e. 
values of the both functions Q, E are the same at these points), then this point is clearly 
a singular point of the curve, namely, a self-crossing (or double) point. This situation 
is illustrated in Fig. 3. In short, we say that the map from the Schottky double to the 
curve glues some points together. This means that E in general can not be realized as a 
smooth algebraic curve in C 2 . 

Let us give a more precise characterization of the points that can be glued together. 
Two distinct points P U P 2 E E are glued together if Q(P 1 ) = Q(P 2 ) and E{P 1 ) = E(P 2 ). 
As is readily seen, this is possible only if Pi and P 2 are on different sides of the double. 
Let Pi = (zi, +), P 2 = (z 2 , — ) with z lj2 G D u. Then these points are glued together if 
S(zi) = z 2 and S(z 2 ) = z\. It might seem that the second equality follows from the first 
one by virtue of the unitarity condition. However, the unitarity condition holds for the 
whole algebraic function S(z) while S(z) in the both equalities above means one and the 
same branch of the Schwarz function obtained as the analytic continuation of z to D oil . 
So, we have two independent conditions to determine the points z± t 2- To put it differently, 
z\ (or z 2 ) is a point in D oi i where two (or more) branches of the Schwarz function take 
same value. 

Nevertheless, the (smooth) Schottky double always has a realization as a smooth 
complex curve P in a higher dimensional space. This curve can be obtained from P as 
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Figure 3: From a smooth Schottky double to a singular curve r and then to the smooth 
curve r. 



a result of the desingularization process. In the case of a self-crossing point this simply 
amounts to incorporating a third meromorphic function on S, say K(P), such that it 
takes different values at the points to be glued together by the map P i— > (Q(P), E(P)) 
(see Fig. 3). The curve is then defined by a system of polynomial equations Ri(Q, K) = 0, 
R 2 (K,E) = in C 3 . If this curve is still singular, the process should be repeated until 
one arrives at a smooth curve t. 

The so obtained smooth curve r is isomorphic to the Riemann surface of the Schwarz 
function. The latter can be also realized as a branched covering of the z-plave. In this 
realization, the antiholomorphic involution is the Schwarz reflection 

z^~S\z) 

which sends a point z to its "mirror" with respect to the boundary contour 3 . 
3.3 The generating differential 

The explicit form of the generating differential for algebraic domains is: 



To prove this, one should check equations (6) for all vector fields d T ( a ) . These equations 
immediately follow from the behaviour of singularities of the Schwarz function under 
Laplacian flows. 

Proposition. The location of poles of the Schwarz function of any algebraic domain 
remains unchanged under Laplacian flows with point-like pumps. The residues remain 

3 This geometric interpretation requires some care. If the point z is not far from the contour, so the 
image belongs to the same sheet of the covering, this is indeed a kind of reflection in the usual sense. 
However, if z is far enough from the contour, the image goes to another sheet across a branch cut. In this 
case the unitarity condition (involutivity of the reflection) is fulfilled only after taking proper branches 
of the Schwarz function. 




S(z)dz on the front side 



zdS(z) on the back side 
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unchanged as well, except for the residue at the pole where the pump is placed: 

d T(b) res z=a (S(z)dz) = l \ | ^ (10) 



Remark. If S(z) is initially regular at a point a, the Laplacian flow with the sink at a 
generates a simple pole at this point with the residue —T^ a \ 

One of the methods to prove this proposition is to consider the function 



C{z) = f 



whose variation under Laplacian flows can be found by a direct calculation. On the other 
hand, as it follows from the properties of Cauchy integrals, singularities of the Schwarz 
function outside 7 (i.e. in D oil ) coincide with singularities of the analytic continuation of 
the C(z) from inside to outside 7. 

It follows from the proposition that the only singularity of the differential d T ( a ) (QdE) 
on the front side is the simple pole at a with residue —1. The unitarity condition means 
that this differential is purely imaginary on 7. Indeed, acting by d% and dr on the 
unitarity condition S(S(z,T), T) = z, we get the identities 



S'(S(z)) S'(z) = 1 (11) 



d T S(S(z), T) + S'(S(z)) d T S(z, T) = (12) 
(here 8t means any of d T ( a )). Combining them, we obtain: 



d T S(z,T)dz = -d T S(S(z),T)dS(z) 



which yields dxSdz = —drSdz on 7. Taking into account the normalization condition, 
we conclude that 

d T(a) (QdE) = -dW {a) for all a e D oil (13) 

Let us fix a and set T = T^, W = . If one takes the function W = W{z,T) as 
a connection in the space of real Riemann surfaces, instead of the function z, i.e. keeps 
W constant when applying the vector field dx, then equation (13) acquires the form 

d T E(W, T)dQ - d T Q(W, T)dE = dW (14) 

In the simply-connected case, it turns, after the restriction to 7 and the identification 
E = f , Q = f, iW = 0, into the Laplacian growth equation (2). 



4 Whitham equations for Laplacian growth of alge- 
braic domains 

In this section we represent the Whitham hierarchy (5) as a finite system of differential 
equations for branch points of the Schwarz function and give a general form of the 
generalized hodograph solution of this system. 
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4.1 Branch points of the Schwarz function 

Given an algebraic domain, consider zeros of the differential dE on the Schottky double. 
Since E = z on the front side, all these zeros are on the back side. Let us denote them by 
XI = (Afc, — ), so Afe's are their projections on the z-plane. Since = dE(\* k ) = dS(Xk), 
Afc's are zeros of the derivative of the Schwarz function in D oil : 

S'(X k ) = , X k e D oil 

Given poles of the Schwarz function and the number of droplets, the number of such 
points can be found from the following reasoning. The differential dE has poles of orders 
rij + 1 at the points (pj, — ) on the back side of the double (recall that pj are poles 
of the Schwarz function) and the second order pole at infinity on the front side. The 
total number of poles is 2 + J2j( n j + 1)) counting multiplicities. For any meromorphic 
differential on a smooth genus-g Riemann surface the difference between the numbers of 
its zeros and poles equals 2g — 2. Therefore, the number M of zeros of dE is given by 

M = 2<7 + 5>,- + l) (15) 

j 

In particular, if all poles pj of the Schwarz function are simple (rij = 1), then 

M = 2(g + N) (16) 

In general position zeros of the differential dE are all simple. Assuming this, the 
critical values of the function E, i.e., values of E at the zeros of dE, 

E k :=E(X* k )=S(Xk) (17) 

play the role of local coordinates in the space of algebraic domains with given N and g 
(see [12]). They are images of the points X k under the Schwarz reflection. Their geometric 
meaning is clarified by expanding the unitarity condition in a vicinity of the points X k 
(or E k ). As a result, we obtain: 

S(z) = S(E k ) + a k ^jz -E k + ... , z^E k (18) 

Therefore, E k s are branch points of the Schwarz function. In general position they are 
of order two. If all poles of the Schwarz function are simple, then all the 2(g + N) branch 
points (see (16)) are inside water droplets. In this case a single-valued branch of the 
Schwarz function can be defined by making g + N cuts between them. 

4.2 Dynamics of the branch points 

Here we derive differential equations for dynamics of the branch points of the Schwarz 
function Ei = Ei(T). 

Let us expand the equality (5) on £ at the points where E fails to be a good local 
parameter, i.e., at zeros A£ of the differential dE. Near these points the local parameters 
are \JE — E k , so the expansions have the form: 

W (a) = w^ a \X* k ) + Cl a) y/E - E k + 0(E - E k ) 
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and similarly for W^ b \ Comparing the singular parts near the branch points, we obtain 
the equalities 

Cf d T(a) E k = Cf d TW E k , cf dW^(X* k ) = cf dW^(Xl) 

for all k. Dividing one by another, we obtain a closed system of differential equations of 
the hydrodynamic type for the branch points E k : 

d T(a) E k = V k {ab) (E 1 , ...,E M ) d T(b) E k , k = 1, . . . , M (19) 
Here the "group velocities" vjf^ are given by: 



In general they are complicated nonlinear functions of E k . The system (19) is diagonal 
with respect to the E k s, so the branch points of the Schwarz function are Riemann 
invariants. 

Admissible Cauchy data for the system (19) are given by the ordinary differential 
equations of the form 

which is equivalent to the "string equation" (14). Similarly to the calculation given above, 
these equations are obtained by expanding (14) near zeros of dE. Note that drE(W,T) 
at any zero A£ of dE is equal to dxE k . 

Solutions of these equations can be constructed by means of the generalized hodograph 
method [14, 12]. In the version of [12], these solutions are obtained by the representation 
of the generating differential in the form 

QdE = dA- T (a) dW {a) - T (b) dW (b) (22) 

where dA is a differential with T-independent singularities: dr(AA(z)) = (A is the 
Laplace operator). For algebraic domains, dA is a meromorphic differential with fixed 
principal parts at all poles. Since Q is a regular function on the back side of E (except 
for infinity), the r.h.s. of (22) must vanish at the points \* k . This requirement leads to 
the "hodograph relations" 

T {a) dW {a \E k ) +T {b) dW {h \E k ) = dA(E k ) 

or, equivalently, 

T (b) +V (<>b) T (a) = u (b) (23) 

for all k — 1, . . . , M, where 



rAb) ( dA(z) 
k ~ \dWW{ 



■ ] ' z=E k 
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They give an implicit solution to the Whitham equations (19) and to the "string equation" 
(21). The Laplacian growth with a sink at the point b corresponds to changing with 
keeping fixed. 

As a rool, the solutions constructed have bifurcation points. In applications to non- 
linear waves, such a point means that the wave front turns over. In the Laplacian growth, 
the bifurcation points correspond to the unphysical singularities (cusps) of the moving 
interface, which are however typical [1, 8] in the limit of zero surface tension. 



5 Example: N = 1 

Consider an example, where the Schwarz function has a simple pole at a point p G D oi i 
and is regular everywhere else in D oi i (N — 1). Here we follow refs. [20]. In the vicinity 
of the pole we have 



S(z) 



p — z 



+ 0(1), 



p, p e D, 



oil 



where the residue is put equal to — //. According to the general argument, S(z) forN = l 
obeys an algebraic equation 

]T A lk S\z)z k = 

l,k=0 

which is quadratic in each variable. Therefore, the Riemann surface of the Schwarz 
function has two sheets over the z-plane. One of them is "physical". It is the sheet 
where the relation S(z) = z holds on the contour. In what follows, we call it the first 
(or upper) sheet. The pole p is on this sheet. The analytic continuation of the Schwarz 
function inside water droplets has another pole at a point which is the value of the 
function E at the simple pole of the function Q at oo on the back side of the Schottky 
double. This pole, q, lies on the second (unphysical) sheet. To summarize, the poles of 
the Schwarz function are as follows: 



S(z) 



// 



p — z 

V 



q 



+ 0(1), 



+ 0(1), 



p on the 1-st (physical) sheet 



q on the 2-nd sheet 



(24) 



The unitarity condition implies that 



S(z) = 



— + 0(z 2 ) , z — > oo on the 1-st sheet 

z 



z 



(25) 



oo on the 2-nd sheet 



Therefore, all singularities of the differential Sdz on the Riemann surface of the Schwarz 
function are two simple poles at p and q on different sheets (with residues —[/,, —v) and 
two second order poles at two infinities (with residues V and Jx). Sum of the residues 
must be zero, so ji + v is a real number: 



H + v = n + v 



(26) 
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Substituting the principal parts of the Schwarz function (24) into the equation 



z 2 z 2 + bz 2 z + bzz 2 + czz + dz 2 + dz 2 + ez + ez + h = (27) 

where z = S(z), we can fix all the unknown coefficients except for the free term h (cf. 
[21]): 

b = -p-q 

c= \p + q\ 2 + pL + v 

d = pq 

e = —p q (p + q) — pv — q/i 

(it is implied that /i, v ^ 0). 

How many droplets may exist for N — 1? Equation (27) in general defines a smooth 
curve of genus 1 (an elliptic curve). The real section of this curve (the set of points such 
that z — z) consists of at most two disjoint closed contours on the plane. If there are 
two contours, a detailed analysis shows (see, e.g., the proof of Theorem 7 in [19] or the 
corresponding example in [20]) that one of them necessarily lies on the unphysical sheet. 
Therefore, two droplets are impossible for N — 1. In other words, all algebraic domains 
with N — 1 are simply-connected. In this case the curve (27) is a rational (genus-0) curve 
with singular points (a smooth genus- 1 curve of the form (27) can not be birationally 
equivalent to the genus-0 Schottky double of a simply-connected domain). Given the 
coefficients b, c, d, e, this condition determines h. 

Solving the quadratic equation (27), we obtain the Schwarz function in a more explicit 
form: 



S{z) = -{P + q) + ~, T7 : 28 

2 2(z-p)(z-q) 

Here P±{z) is a polynomial of 4-th degree with the leading term z 4 . Generally speaking, 



the function y P^z) has four branch points which are simple roots of the polynomial P4. 
The degeneration of the curve occurs when two roots, say E 3 and £4, merge: 



y/P&j =(z- E 3 )JK(zj, P 2 (z) = (z- E^z - E 2 ) 

so there is only one cut between Ei, E 2 and the curve is of genus 0. The second cut 
degenerates into the double point of the curve. Comparing with (25), one can see that 



the 1-st (physical) sheet corresponds to the branch of y P 2 (z) such that y P 2 (z) — > z as 
z — > 00. 

Endpoints of the cut E 1 , E 2 in general can not be found explicitly in a simple form (in 
the case at hand one has to solve the equation of 4-th degree). We know, however, that 
they obey the equations of the hydrodynamic type and the solutions can be represented 
in the hodograph form (22). 

The hodograph relation follows from the analytic properties of the differential Sdz. 
It is a meromorphic differential on the 2-sheeted Riemann surface with a branch cut 
between E 1: E 2 . It has two simple poles and two poles of second order. The simple poles 
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Figure 4: The Riemann surface of the Schwarz function for N = 1 is a two-sheeted 
covering of the z-plane. The map E : E ^ P = C Li {00} is a two-sheeted covering, i.e., 
each point z G P has two pre-images (counting multiplicities). The image of the front 
side of the double is the domain D u on the first sheet. The image of the back side of the 
double is the rest part of the two-sheeted Riemann surface. Namely, it is the interior of 
the boundary contour on the hrst sheet (the domain D ) and the whole second sheet. 



are p on the 1-st sheet and q on the 2-nd sheet with residues — /i, —v respectively. The 
second order poles are 00^ (00 on the 1-st sheet) and 00^ (00 on the 2-nd sheet) with 
residues V and p, respectively. Therefore, we can write 

Sdz = qdW + +pdW_ - iidW pWj0oW + JidW^^v - {fx- l^dW^,^) (29) 

Here dW± are meromorphic differentials with the only pole of second order at infinity 
on the first and second sheets respectively with the principal part (1 + 0(z~ 2 ))dz as 
z — > 00 and dW a ,b is the dipole meromorphic differential with simple poles at a and b 
with residues ±1 respectively. The upper indices indicate the sheet where the points lie. 
For example, dW ^1)^2) is the dipole differential with poles at infinity on the first sheet 
and q on the second sheet. Explicitly, these differentials are: 



dW± = 



Ei + E 2 



dz 



dW, 



^(z) + jFMa) v /^) + v /^(6)\ dz 



<i.i 



z — a z — b 



2JP 2 (z) 



Plugging the explicit formulas for the differentials into (29), we represent S(z) in the 
form 

where 

2f 1 (z)=p + q + -^— + 



p — z q — z 

2f 2 (z) = -(q-p)(2z-E 1 -E 2 )+p-u + f ,^^ + u^—^ 
2 z — p z — q 

The hodograph equations follow from the fact that the values S(E 1 ) and S(E 2 ) are finite, 
i.e., f 2 (E 1 ) = f 2 (E 2 ) = 0. 
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If all the parameters are real and q < E 1 , E 2 < p, the system of hodograph equations 
acquires the form 



, , , , (30) 

where square roots of positive numbers are assumed to be positive. The real parameters 
[/,, T, where 

™ 1 f „, s , Area(D) 

are local coordinates in the two-dimensional variety of contours obtained from a circle 
(of radius ro) by Laplacian growth processes with oil pumps located at p and oo. In this 
interpretation, T — 7q is the total amount of sucked oil and p is the amount of oil sucked 
from the point p. 
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